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§ ■ ABSTRACT 

O 
(N 

j>^. Recent observations of M87 and some blazars reveal violent activity in small 

regions located at relatively large distances from the central engine. Motivated 
by these considerations, we study the hydrodynamic collimation of a relativistic 
cooling outflow using a semi-analytical model developed earlier. We first demon- 
^ . strate that radiative cooling of the shocked outflow layer can lead to a focusing of 

Q^. the outflow and its reconfinement in a region having a very small cross-sectional 

Q . radius. Such a configuration can produce rapid variability at large distances from 

. the central engine via reflections of the converging recollimation shock. Fossi- 

ck . ble applications of this model to TeV blazars are discussed. We then apply our 

model to M87. The low radiative efficiency of the M87 jet renders focusing un- 
^ , likely. However, the shallow profile of the ambient medium pressure inferred from 

^ . observations, results in extremely good collimation that can explain the reported 

if^ '. variability of the X-ray flux emitted from the HST-1 knot, 

p , 
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j_j ■ 1. Introduction 

The interaction of a relativistic jet with the surrounding matter is likely to play an 
important role in its collimation and the dissipation of its bulk energy. Such an interaction 
is evident on the largest scales, where hot spots, presumably associated with termination 
shocks, and features that appear to be associated with recollimation by the surrounding 
matter and by backflows from the jet endpoints are often observed in AGNs and, in some 
cases, microquasars. Less clear is the effect of the environment on relativistic jets at much 
smaller scales, in particular in the blazar zone. 
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Recent observations of TeV AGNs suggest that a considerable fraction of the bulk en- 
ergy may dissipate in reconfinement shocks on VLBI scales. These observations motivate 
reconsideration of the standard view, according to which the broad-band, highly variable 
emission seen in blazars is produced predominantly behind internal shocks that form in 
colliding fluid shells. One indication is the systematic differences between the large values 
of the Doppler factor inferred from TeV observations (e.g., Levinson 2006; Begelman et 
al. 2008) and the much lower values inferred from unification models (Urry and Padovani 
1991; Hardcastle et al. 2003) and superluminal motions on parsec scales (e.g., Marscher 
1999; Jorstad et al. 2001). Such differences can be naturally accounted for if dissipation 
occurs in quasi stationary patterns, e.g., oblique shocks, that result from the interaction of 
the jet with ambient matter. In this scenario radio observations reflect the pattern speed 
while the Lorentz factor inferred from TeV observations is associated with the speed of the 
fluid passing the structure. Furthermore, the large Doppler factors implied by opacity argu- 
ments may not be required in the first place if the TeV emission is produced at radii larger 
than commonly believed. However, in that case channeling of the bulk energy into a rather 
small area needs to be assumed in order to account for the rapid variability often observed. 
Alternative explanations have also been offered in order to resolve the so called 'Doppler 
factor crises', including a structure consisting of interacting spine and sheath (Ghisellini et 
al. 2005), opening angle effects (Gopal-Krishna 2004) and jet deceleration (Georganopoulos 
& Kazanaz 2003; Levinson 2007). 

Based on analysis of multi-wavelength observations during the 2005 outburst in 3C454.3, 
Sikora et al. (2008) argued that the blazar emission in this source is produced at a distance 
of several parsecs from the putative black hole, where the mm photosphere is located. They 
proposed that the blazar activity is driven by a standing reconfinement shock that, they 
claim, dissipates energy more efficiently than internal shocks. They further argued that the 
high energy emission is likely produced via inverse Compton scattering of IR dust emission, 
and demonstrated that the observed SED can be reconstructed in such a model. A fit of their 
model to the observed broadband spectrum favors a relatively large Lorentz factor of the 
emitting fluid, F ~ 20, consistent with that required to account for the temporal variations 
of the optical and millimeter fluxes. 

Perhaps the best example that an outflow energy can be channeled into a small area 
located far from the central engine is the HST-1 knot in M87, a stationary radio feature 
associated with the sub-kpc scale jet. The knot is located at a projected distance of 60 pc 
(0.86") from the central engine, and is known to be a region of violent activity. Sub-features 
moving away from the main knot of the HST-1 complex at superluminal speeds have been 
detected recently (Cheung et al. 2007). In addition, rapid, large amplitude variations of the 
resolved X-ray emission from HST-1 have been reported, with a doubling time iyar of 0.14 yrs. 



-3- 



The observed variability limits the linear size of the X-ray source to d ^ Soctvaj: ~ 0.04A5d pc, 
where is the Doppler factor, which for a viewing angle of 6n ~ 30° is at least three orders 
of magnitude smaller than the distance between the HST-1 knot and the central black hole. 
The constraints on the size of the radiation source are much more severe if the observed TeV 
emission also originated from HST-1, as suggested by some based on a claimed correlation 
between the X-ray and TeV emission (e.g., Cheung et al. 2007, but c.f., Neronov & Aharonian 
2007). It has been proposed that HST-1 is associated with a recoUimation nozzle (Stawarz 
et al. 2006; hereafter ST06). In this picture the rapid variability and superluminal sub-knots 
that seem to be expelled from the HST-1 complex can be associated with shocks produced 
by reflection of the recoUimation shock at the axis (Levinson and Bromberg 2008). The 
rapid variability sets a limit on the cross-sectional radius of the channel at the location of 
HST-1 that depends on the fraction of jet power radiated as X-rays (and TeV emission, if 
indeed originating from the same location). The conditions required to produced a structure 
consistent with this picture are examined in §3 below. 

In what follows we exploit a model developed by us earlier to address some of the issues 
discussed above. We generalize it to include radiative losses and demonstrate that even 
modest radiative cooling of the shocked jet material can lead to extremely good focusing. 
The effect of cooling on the coUimation and confinement of non-relativistic jets has been 
studied earlier in the context of SS433 (e.g., Peter & Eichler 1995). Our work presents 
an extension of this idea into the relativistic regime. A preliminary account of the results 
presented below is outhned in Levinson and Bromberg (2008). 



2. The Model 

In a previous paper (Bromberg and Levinson 2007; hereafter BL07) we constructed a 
class of semi- analytical models for the confinement and coUimation of a relativistic jet by 
the pressure and inertia of a surrounding medium. Both, confinement by kinetic pressure of 
a static corona, and confinement by the ram pressure of a supersonic wind emanating from 
a disk surrounding the inner source have been considered. In general, the collision of the 
inner jet with the confining medium leads to the formation of a contact discontinuity across 
which the total pressure is continuous, and an oblique shock across which the streamlines of 
the colliding flow are deflected. In cases where conflnement of the inner flow is accomplished 
through collision with a supersonic wind a second shock forms in the exterior wind. The 
model outlined in BL07 computes the structure of the shocked layers of the deflected inner 
jet and the exterior wind in the latter case, assuming a steady, axisymmetric flow. In BL07 
the focus was on the application to GRBs. Radiative losses have been ignored since the 
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large optical depth of the shocked jet layer on scales of interest renders such losses negligibly 
small. In blazars, recollimation shocks are expected to form above the photosphere. If the 
cooling rate of the shocked plasma is high enough, then a significant fraction of the thermal 
energy may be radiated away and this may affect the structure of the flow. To study this 
effect we incorporated radiative losses into our model. The results are described below. 



2.1. Analytical approximations 

Consider the interaction of a relativistic outflow with a gaseous condensation extending 
from some fiducial height Zq above the equatorial plan to infinity, and having a pressure 
profile Pext{z) = Po{z / Zo)~^ . The outflow is assumed to be ejected from a point source into 
a cone of opening angle 6j, with a total power Lj distributed uniformly inside the cone, 
and velocity Pjo at z = zq. As explained above, the structure of the confined flow consists 
of a contact discontinuity separating the shocked jet layer and the ambient medium and 
a collimation shock 0. The details of this structure will depend, quite generally, on the 
parameters of the injected outflow and on the external pressure. We denote by rc{z) and 
rs{z) the cross-sectional radii of the contact discontinuity and shock surfaces, respectively. 
By employing the shock jump conditions and requiring momentum balance at the contact 
discontinuity surface, the hydrodynamic equations can be reduced to a coupled, nonlinear 
set of ODEs for rc{z), rs{z) and the thermodynamic parameters (density, temperature, etc.) 
of the shocked fluid (BL07). In general, solutions must be sought numerically. However, the 
problem can be considerably simplified in certain regimes. At small angles one obtains to 
lowest order (Komissarov & Falle 1997; BL07), 

±l = !l-Atane,z'-'^/\ (1) 
dz z 

the solution of which is 

r,{z) = ztanOj - ^^ztanOj [z^-'^l'' - zj"''/') , (2) 

where A = {ttcpoZq / Lj^jo^iY^"^ . Here = {1 — n ■ f3j^/h ■ (3-) ^ 1, with /3_ and (3+ denoting 
the local fluid 3-velocity upstream and downstream of the shock, respectively, and n the 
normal to the shock surface. BL07 adopted = 1 for convenience whereas ST06 adopted 
^1 = 0.7. The point z* at which the shock reaches the axis is determined from the condition 



^For certain external pressure profiles the flow equations admit self-similar solutions with no shocks 
(Zakamska et al., 2008). 
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rs{z = z*) = 0. Using eq. ([2]) one finds, 

^-(^^^» ) ■ 

Note that z* depends solely on the total jet power Lj and the external pressure profile. 
However, as will be shown below, the jet profile depends also on the opening angle and other 
details. From eq. ([3]) it is readily seen that convergence is guaranteed for < 2. In the case 
T] > 2 the shock will approach the axis provided po exceeds some critical value pc, given by 

. (VzfhiE. ,4) 

Ancz^ 

As Po approaches pc the location of the refiection point z* approaches infinity, and when 
Po < Pc the shock diverges. The structure of the polar outfiow in the latter case consists 
of a core containing the unshocked jet enveloped by the shocked jet layer that expands 
relativistically, but with a Lorentz factor considerably smaller than that of the unshocked 
jet, as described in detail in BLOT. 

Approximate solutions can also be obtained analytically for the contact discontinuity 
surface. In the absence of radiative energy losses the radius rc{z) at 2; >> is given 
approximately by (BLOT) 

r,(z) = rcoiz/zor^\ (5) 

Evidently, collimation occurs when rj < A, however, the fiow never shrinks, even in situations 
where the shock itself reaches the axis. The reason is that the pressure inside the shocked 
layer pushes against the ambient gas and keeps it expanding as long as > 0. The situation 
may change if the shocked gas cools radiatively at a high enough rate that allows strong 
compression shocked jet material. In that case, the shocked layer remains thin, leading to 
the convergence of the contact discontinuity towards the axis. In the limit of extremely 
rapid cooling the width of the shocked layer approaches zero, viz., Tc — Vg, and focusing is 
expected under the same conditions that lead to convergence of the collimation shock. 



2.2. Radiative reconfinement and focusing of a relativistic flow 

A fraction of the shocked energy is likely to be tapped for acceleration of particles 
to nonthermal energies, e.g., by shock acceleration, or shear mechanism (Ostrowski, 1998). 
The thermal and nonthermal electrons may cool via synchrotron and/or inverse Compton 
emission, leading to compression of the shocked layer. This compression is clearly inde- 
pendent of the details of the cooling mechanism, but merely on the fraction of the total 
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Fig. 1. — Lorentz factor of the shocked and unshocked jet in the case of confinement by a static 
corona with a pressure profile pext oc z"^, for different values of the cooling parameters. In all cases 
shown = 0.01. The ratio of the total luminosity radiated away by the shocked jet layer and 
total jet power is Lc/Lj = (no radiative losses) in the upper panel, Lc/Lj =0.1 (^e = 0.015) in 
the middle panel, and Lc/Lj = 0.27 (.^e = 0.06) in the lower panel. The injected flow in all panels 
consists of a cold, purely baryonic fluid with Lorentz factor Tjo = 10 at the injection point z = zq. 
The shock surface is marked by the dashed line 

kinetic energy (the sum of thermal energy and the energy carried by the nonthermal pop- 
ulation) that is being lost to radiation over dynamical scales. We suppose that a fraction 
C,B of the total energy flux in the shocked layer is carried by magnetic fields. The pa- 
rameter can be readily related to the sigma parameter of the injected outflow via the 
shock jump conditions. We consider sufficiently low sigma flows for which the magnetic 
field is dynamically unimportant. Let us assume first that the shocked plasma is in a ther- 
mal state, and that electrons and protons are strongly coupled (that is, the equilibration 
time is much shorter than the flow time). If the specific enthalpy in the shocked layer 
is dominated by rest mass, as we find to be the case for cold, baryonic jets, then to a 
good approximation U'^ = B''^/8n = ^B^pC^^p (the prime refers to quantities measured 
in the rest frame of the flow). Since the shocked layer is confined by the external pres- 
sure the average energy of thermal electrons is e'^ = 'j'^meC^ = pextf^n'^- With n'^ = n'^ 
and Pext = Po{z / zo)~^ , the ratio of cooling time due to synchrotron and inverse Compton 
emission, t'^ooi = 5 x 10^i?'~^7g~"'^(l + UL/U'b)~^ s, where UU is the energy density of the 
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Fig. 2. — The dependence of the solution on jet temperature. The external pressure profile and the 
value of S^B are the same as in fig. 1. The parameter was adjusted such that Lc/Lj = 0.27 in all 
panels. The power of the injected flow is given to a good approximation by Lj = h'MjTcP, where 
h' denotes the proper entropy per baryon of the unshocked fluid. The ratio of total jet energy to 
rest mass energy is Lj/MjTc^ = h' = 10 in the upper panel, Lj/MjTc^ = 4 in the middle panel, 
and Lj/MjTc^ = 1 in the lower panel. 



background radiation field measured in the fluid rest frame, and flow time, = z/cT, is 
given by: 

^_U = 2(i + ^V (^] " ( ^ 1 " (^] " (^) r (6) 

t'j \ U'b) VO-ly dyn cm-2y VO-1 pc; \zo ) ' ^' 

Thus, cooling of the thermal plasma on timescale shorter than the flow time requires Up^/U'^ 
~ (^B/0.1)~^r for the above choice of parameters. If a fraction of the disk luminosity 
Ld is scattered across the jet by the surrounding gas (e.g., Blandford and Levinson 1995), 
then Uph — 'TdLd/4:7!'cr'^ , and t/^^^ = T^Uph on account of beaming of the radiation field in 
the comoving frame. Rapid cooling of thermal electrons then implies L^/Lj ^ O.lr^^F^^^^^, 
where 6j is the opening angle of the jet, and may be expected in ERG blazars. 

It could well be that a small fraction of the electrons in the shocked layer are being 
continuously picked up from the thermal pool and injected to nonthermal energies. This 
population cools rapidly and can significantly enhance the radiative loss of the shocked layer. 
For instance, over- stability of the contact discontinuity surface (e.g., to Kelvin-Helmholtz 



- 8 - 



r 

4 5 6 7 8 9 




z/r 

CO 



Fig. 3. — The effect of the ambient medium. In aU cases shown FjQ = 10, Lj/TMjC^ = 1, = 0.01, 
and = 0.04. The upper, middle and lower panels correspond to r/ = 1, 2, and 3, respectively. 
The radiative efficiency L^/Lj is indicated in each panel. 

modes) may lead to nonlinear oscillations of the contact surface and the reconfinement 
shock itself that can enhance local dissipation in the region between the contact and the 
reconfinement shock (e.g., via generation of steepening waves, or by locally changing the 
incidence angle of the fluid upstream of the reconfinement shock), giving rise to acceleration 
of electrons and positrons to nonthermal energies. Or it could be that a fraction of the 
particles are Fermi accelerated in the sheared flow as described in Ostrowski (1998). In order 
to account for such processes we suppose that a fraction of the thermal energy behind the 
reconfinement shock is injected as a power law distribution of electron^: dn'^/de'^^ oc e'~'^. 
For simplicity we consider only synchrotron cooling of the nonthermal population. Fixing 
and we then compute the synchrotron emissivity of the shocked gas, and incorporate 
it into the flow equations of the shocked layer. The details are given in the appendix. To 
simplify the analysis further the fractions and ^b are taken to be constants along the 
channel (independent of z). The free parameters of the model are the power Lj, opening 
angle 9j, Lorentz factor Fjo = '^j{z = zq), the mass flux Mj and the product ^e^^'^ the 
injected flow, and po, rj of the external medium. 



^Tlie sliape of the electron spectrum depends on the details of the acceleration mechanism at work. Since 
we are merely interested in the net cooling rate the use of a single power law is sufficient. 
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Examples are shown in figures 1-3. In all the examples shown the confining medium 
is static with a pressure profile Pext = Po{z/zq)~^, as discussed above, and the jet power is 
Lj = 10"^^ erg/s. The integration starts at 2:0 = 1 pc, where the opening angle of the injected 
flow is taken to be 6' = 10°. Integration of the full set of equations confirms the scaling 
derived analytically in eq. Q. The effect of cooling is examined in fig. 1. In this example 
the ejected flow is assumed to be cold, in the sense that its energy density is dominated by 
rest mass energy. To be more precise, we take Lj = MjTjQC^ at the initial impact point 
z = zq, (as stated above, in all cases considered here << 1 and has a little effect on 
the shock jump conditions that we ignore). We considered both, baryonic jets and pure 
electron-positron jets and, as expected, found little dependence on the outflow composition, 
with the exception of the temperature behind the shock which is higher in the baryonic case, 
owing to the lower density. The difference between the three cases shown in flg. 1 is in the 
synchrotron cooling rate behind the oblique shock, which is controlled by the parameters 
and ^B- The ratio of the total luminosity radiated away by the shocked jet layer. 



where is given in Eq. (fT2l) in the appendix, and total jet power Lj in the upper, middle 
and lower panels is Lj = (no radiative losses), 0.1, and 0.27, respectively. As seen, in all 
cases the shock (indicated by the dashed line) approaches the axis at the same distance from 
the injection point, consistent with the scaling in eq. Q. The main effect of the radiative 
cooling, as illustrated in flg. 1, is to increase the shock compression ratio, thereby reducing 
the width of the shocked layer and, as a result, the cross- sectional radius of the jet at the 
point z*. As demonstrated in the bottom panel, substantial focusing can be accomplished 
even for a modest radiative efficiency. We have also made some runs where only cooling of the 
thermal plasma has been incorporate into the model (^e = 0) and found the same behavior 
essentially. Signiflcant focusing occurs, as expected, when tcooi ~ in Eq. In reality, 
the level of focusing may be limited by other components that have been neglected in our 
model, e.g., magnetic pressure (see further discussion below), or the pressure of nonthermal 
baryons accelerated by the sheared flow. 

The effect of cooling on the structure of the flow is less dramatic in hot jets, in which 
the dimensionless enthalpy per particle of the unshocked fluid is much larger than unity (see 
flg. 2). The reason is that in that case the bulk energy of the shocked fluid is dominated by 
pressure rather than rest mass energy and, therefore, a larger fraction of the total jet energy 
must be radiated away in order to signiflcantly compress the shocked layer. 

Figure 3 exhibits solutions obtained for different proflles of the external pressure. In 
this example we flxed the parameters ^b, ^e, Lj, Mp and po and varied rj. Fixing ^b and 




(7) 
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Lj means essentially that the magnetic field in the shocked layer is the same for all cases. 
As expected from eq. ([3]), the location of the refiection point is closer for smaller values of 
rj. This is a consequence of the larger momentum transfer across the interface separating 
the ambient medium and the jet, that forces larger defiection angles of streamlines crossing 
the shock. Since the transverse component of the incident energy fiux dissipates behind the 
shock this also means larger dissipation. This is the reason why the Lorentz factor of the 
shocked fiuid is smaller for smaller values of t], as seen in fig. 3. The larger dissipation 
also leads to a much faster cooling and, as a result, a much better focusing of the jet. The 
fraction of the jet power that has been radiated away in each case is indicated. 

The above analysis neglects the contribution of magnetic pressure in the shocked jet 
layer. This pressure may support the cooling layer and suppress focusing in cases where 
the field is predominantly poloidal or turbulent, and strong enough. On the other hand, a 
toroidal component gives rise to magnetic hoop stress that actually helps coUimating the 
fiow. The presence of a helical component has been inferred from radio polarization maps 
in some blazars (e.g., Gabuzda et al. 2004). At any rate, we have verified that in all cases 
studied above the magnetic pressure was not important, except near the point z* in the 
upper panel in Fig 3. 

In steady state, reflections of the converging shock at z* (where the shock crosses the 
axis) generally leads to formation of conical shocks. However temporal fluctuations may lead 
to more complicated structures, in particular formation of "internal" shocks in the vicinity 
of the recoUimation nozzle. To model such effects requires 2D, time dependent simulations. 
There is evidence for such features in the full 2D simulations performed by Alloy et al. 
(2005). A considerable fraction of the remaining jet power may dissipate via these internal 
shocks in a region much smaller than what would be expected in the case of a conical jet. 
This may lead to large amplitude variations over timescales much shorter than the jet radius, 
as discussed further below. 

An important question is whether the narrow structure formed beneath the reflection 
point is stable. Our model cannot account for any temporal effects, particularly those asso- 
ciated with Kelvin-Helmholtz instabilities at the interface separating the shocked jet layer 
and the ambient matter. If the instability grows to a nonlinear state over the expansion 
time then pinching of the fast jet and mixing of ambient material with the jet fluid near 
the interface is anticipated. This can modify the transverse structure of the jet. Numerical 
simulations exhibit some evidence for such instabilities (e.g., Aloy, 2005; Hardee & Hughes 
2003), but the details should depend on the conflguration of the conflning medium, and on 
the structure of the viscous boundary layer. Detailed analysis of the stability of the cooling 
jets considered here requires 3D simulations with an appropriate setup. 
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3. The case of HST-1 in M87 

The rapid variability of the resolved X-ray emission from the HST-1 complex sets a 
limit on the cross-sectional radius of the jet at the location of HST-1 that depends on the 
fraction of jet power radiated as X-rays (and TeV emission, if indeed originating from the 
same location). Estimates of the jet power in M87 yield Lj ^ 10^'^ erg s^^ (Bicknell & 
Begelman 1996). For the observed X-ray (Cheung, et al. 2007) and TeV (Aharonian, et al. 
2006) luminosities, Lxev ~ i^x ~ 10"^^ erg s~^, this implies a rather small radiative efficiency, 
Cr = L^/Lj ~ 10^^. In order to account for the observed luminosity the size scale of the 

1 /2 

fluctuations producing the X-ray flare must satisfy d ^ 2er a, where a denotes the cross- 
sectional radius of the jet at the dissipation region. The variability time, on the other hand, 
constrain the size of these fluctuations to be < ^Dct^ar- Combining the two constraints 

— 1/2 

yields a < 0.5er Soct^sLr- The apparent speed measured for the superluminal sub-knots in 
the HST-1 complex, Papp ~ 4, implies F > 4 for the fluid passing the HST-1 knot and a 
viewing angle 9n ~ 30° . Adopting 6'„ = 30°, F = 4, and = 10"'^, we estimate a < 1 pc for 
the reported X-ray variability, consistent with the associated HST source, and a < 0.05 pc 
for the TeV variability if indeed associate with HST-1 (Cheung et al. 2007). For a conical 
jet this requires an opening angle 9j = a/zusTi < 10^^ and < 5 x 10"^, respectively, that 
seem somewhat unrealistic. 

ST06 proposed that HST-1 reflects the location of the point where the reconflnement 
shock reaches the axis. The picture they envisaged is that, at a distance zq of several parsecs 
from the central engine the jet encounters a gaseous condensation having a pressure proflle of 
the form Pext{z) = Pb{z/zb)~^, with the values zb ~ 200 pc, pb = 1-5 x 10"^ dyn cm~^ and 
rj = 0.6 adopted in the reconflnement zone. By employing a model developed by Komissarov 
and Falle (1997), they computed the proflle of the converging reconflnement shock and found 
that it reaches the axis at a distance 



4 hcpoz'b 



1/(2-^) 



(8) 



from the central engine, which for t] < 2 coincides with eq. ([3]) in the limit Zq 0. For the 
parameters adopted by ST06 a total jet power of Lj ~ 10"^^ ergs s~^ is required to match z* 
with the location of HST-1, consistent with the estimate of Bicknell & Begelman (1996) [§. 
The treatment of ST06 is incomplete, as it does not compute the structure of the shocked 



■^STOG imposed the additional restrictions that the kinetic pressure of the jet at z = zq equals the ambient 
pressure, viz., pj{z = zq) = Pext{z = zq), and that the opening angle of the jet at z = zq equals the inverse 
of its Mach number. This uniquely determines Lj and z*. Other situations can be envisaged that do not 
require such restrictions. 
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Fig. 4. — Lorentz factor of the shocked and unshocked layers of a baryonic jet in the case of 
confinement by a static condensation with a pressure profHe pextiz) = pb{z/ zb)~^ ■ The integration 
starts at = 4 pc. In ah cases shown Lj = 5 x 10^^ erg s~^, Tjo = Tj{z = zq) = 10, zb ~ 200 
pc, pb = 1-5 X 10^^ dyn cm~^ and rj = 0.6. The semi-opening angle of the injected outflow is: 
9j = 0.05 (upper panel) 9j = 0.1 (middle panel), and 9j = 0.2 (bottom panel). 

jet layer, and in particular the cross-sectional radius of the jet at the location of HST-1. To 
study this we applied our semi- analytic model to M87. Since only a small fraction of the 
total jet power is converted to radiation in this object, viz., Lc/Lj < 10~^, it is not expected 
to affect the dynamics of the system (see fig. 1). We therefore set C,e = 0. We emphasis that 
even though unlikely to affect the structure, synchrotron emission may still be observable 
from regions where the recoUimation shock is sufficiently strong. 

An example is shown in fig. 4, where the above values of the ambient gas parameters 
have been invoked. In this example the integration starts at 2:0 = 4 pc, where 6j is fixed. 
Starting the integration at smaller distances Zq does not change the solution much. This is 
expected since the ratio of the outflow ram pressure and the external pressure increases as z 
decreases (roughly as z~^''^), so that a.t z « z* the shocked jet layer is very thin and the jet 
profile is essentially unaffected by the ambient matter. As seen the location of the reflection 
point z* is the same in all cases but the profile of the contact discontinuity depends on the 
opening angle of the jet at the initial location zo- It is also seen that coUimation occurs on 
intermediate scales (20 to 40 pc for the parameters used in fig. 4) where the recoUimation 
shock is strongest. Observable synchrotron emission may be produced behind the shock in 
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this region. Beyond the coUimation radius the shock becomes much weaker and we expect 
httle or no radiation there. This appears to be consistent with morphology of the radio jet 
on sub-kpc scales. 

4. Discussion 

We have constructed a class of semi-analytical models for the coUimation of a relativistic 
cooling jet by the pressure of a surrounding medium. To illustrate the effect of cooling we 
assumed, for computational convenience, that radiative losses behind the coUimation shock 
are dominated by synchrotron emission of accelerated electrons. However, the structure of 
the shocked jet layer depends merely on the net cooling rate and not on the details of the 
cooling mechanism, provided momentum losses are small. The main conclusion is that, under 
certain conditions, radiative cooling of the shocked jet layer can lead to a good focusing of 
the jet. Focusing requires ambient medium with sufficiently flat pressure profile, and is most 
effective when the jet is cold, that is, when its energy flux is carried predominantly by rest 
mass. The reason is that in the latter case modest radiative efficiency is already sufficient 
to reduce the temperature behind the coUimation shock considerably, giving rise to a high 
compression of the shocked jet material. The location where the reconfinement shock reaches 
the axis depends solely on the external pressure and the jet power. However, the profile of 
the jet and, in particular, its cross-sectional radius at the reflection point depend also on 
its opening angle and speciflc enthalpy at the base of the reconflnement zone. Our analysis 
does not account for the effect of radiative drag, that might be important in cases where 
ERG emission dominates. The latter may somewhat alter the dynamics and structure of the 
shocked layer, but should still lead to focusing under the same conditions. In fact, we have 
shown that if the intensity of the background radiation field is as large as typically assumed 
in ERG models than cooling of the thermal electrons alone may be rapid enough to cause 
substantial focusing. Since the Lorentz factor of the flow near the contact discontinuity 
surface is smaller than that of the unshocked jet the emission from this region is expected to 
be less beamed. The GeV emission detected recently by Fermi in the radio galaxies Gen A 
and NGG 1275 (Abdo, et al. 2009a,b) may originate from this region. Temporal fluctuations 
of the central engine may lead to a more complex behavior. In particular, internal shocks 
may form below the point where the coUimation shock intersect the axis and contribute to 
the emission. There is some evidence for such a behavior in BL Lac (Marscher et al. 2008). 

Radiative reconfinement can channel the outflow into a region having a very small cross- 
sectional radius. Reflections of the converging recoUimation shock there should give rise to 
additional dissipation via formation of internal shocks. This can naturally explain ejections 
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of superluminal radio sub-features and rapid variability from an otherwise quasi- stationary 
region, as occasionally observed in radio jets of blazars (Jorstad et al. 2001) and, most 
notably, in M87. The pattern speed of the main radio knot may be associated with the 
location of the reflection point while the Lorentz factor inferred from the variability of the 
VHE emission is associated with the speed of the fluid passing that point. The superluminal 
sub-features ejected from the main knot can be interpreted as reflections of the internal 
shocks that formed at the reflection point. 

The confining medium in blazars may be associated with matter in the broad line region. 
The characteristic pressure measured in broad line emitters is pblr ~ 10~^ — 10~^ dyn cm~^, 
on scales zblr ~ 0.01 — 1 pc. The nature of the broad line emitters is yet an open issue. 
If consists of small clouds, as envisioned by some, than the question remains as to how 
these clouds are confined (if at all). Most likely, some intercloud medium with a similar 
pressure is present. Choosing for illustration po = Pblr, = zblr and 77 = 1 in eq. ([3]), we 
obtain: z^ ~ 2.5(Lj/10'^^erg s~^)(pBiij/10~^dyn cmr'^)^^{zBLR/^-'\- pc)^^ pc. The pressure 
distribution in the BLR is uncertain. However, as illustrated in fig. 3 the location of the 
intersection point, z*, does not depend strongly on rj. Thus, we conclude that reconfinement, 
and even substantial focusing if radiation losses are significant, can occur on VLBI scales, 
as, e.g., proposed for 3C454.3 by Sikora et al. (2008). For the above numbers we find that if 
about 30 percents of the bulk energy is radiated away behind the shock, viz., Lc/Lj ~ 0.3, 
then the cross-sectional radius of the jet at 2* is a ~ 10~^'^2;*. 

In M87 the radiative efficiency appears to be very low and cooling is unlikely to affect 
the structure of the confined jet. However, the fiat pressure profile of the ambient gas on sub- 
kpc scales (ST06) leads nonetheless to extremely good collimation. For the ambient pressure 
profile adopted by ST06 we find a cross- sectional radius of a/rusT ^ 10^^ if the opening 
angle of the jet at the base of the reconfinement zone (~ several parsecs) satisfies 6j ^0.1. 
This level of collimation is sufficient to account for the X-ray variability observed, but not for 
the variability of the TeV emission. The latter would require extremely good collimation on 
sub-parsec scales by some other mechanism. More likely, the TeV emission originated from a 
different location, e.g., the black hole magnetosphere (Levinson 2000; Neronov & Aharonian 
2007). 

Appendix A: Inclusion of synchrotron cooling 

Let Tj'g denotes the stress-energy tensor of the shocked jet layer (including the nonther- 
mal particles). The dynamics of the shocked jet layer is governed by the continuity equation 
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and by energy-momentum conservation: 

a.T/; = 5^ (9) 

where S'^ denotes the source term associated with energy and momentum losses of the 
shocked fluid. Our treatment assumes that energy losses are dominated by synchrotron 
emission of electrons having a power law energy distribution, n'{e'^) = Ke'^"^, with a lower 
cutoff eg which we take to be mgC^ + kTjg for convenience, where Tjg is the temperature 
of the shocked layer, and an upper cutoff 



/2 2 I " ""e"-- / ' acc \ /in\ 

^e.max = + -—^ [^), (10) 



_ ,2 3 m^c^ /ria 

e,max -e- + 2 \ B' 

which was obtained by equating the synchrotron cooling rate with the acceleration rate, 

^acc ~ Vacc^B'/nieC, where r]acc represents the acceleration efficiency, and B' the magnetic 

field, as measured in the fiuid rest frame. The power per unit volume emitted in the fiuid 

rest frame can be readily computed using the synchrotron emissivity jgyni^)- terms of the 

g/ 

total energy density of nonthermal electrons, u'^ — J,'='"""" e'^n'{e'^)de'^, and the magnetic field 
energy density, u'^ — B^/Stt, it can be expressed as: 



1/0 



where T{x) is a Gamma function with argument x. As a closure condition for the set of 
hydrodynamic equations wc take u'^ to be a fixed fraction of the total energy density 
behind the shock and u'^ to be a fraction of the thermal energy behind the shock (which 
for a strong shock is roughly the energy dissipated in the shock). Spcciflcally, Ub = S,BWjs, 
Ue = ie{wjs — njsmc^), where rijs and Wjs are the proper baryon density and enthalpy of 
the shocked gas (see below). We also adopt rjacc ~ 1 for convenience. Since the emission is 
isotropic in the rest frame of the fiuid there are no momentum losses there, and so 5"^ = 
(5"°, 0, 0, 0). Transforming to the star frame we then obtain 

5/^ = AM5'- = r,,5'°(l,/5,-,) (12) 

Here cf3js is the 3-velocity of the shocked fluid and Fjg the corresponding Lorentz factor. 
Using the above derivation, eqs (22) of BL07 is generalized to: 

^ [w,,r,MUrl - r])] - S'%.{rl - r|) = 2^,r,t/,r,£^, (13) 

and eq. (24) to: 
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with S'^ given by eq. f|TT]) . The various quantities are defined in BL07. The total energy 
density of the shocked layer is the sum of the the contributions of the thermal plasma, the 
nonthermal (accelerated) particles, rest mass and magnetic fields, viz., Ujs = fijsmpC^ + 
Ugas + Ue + ub- (In all cases studied above the magnetic field is dynamically unimportant, 
viz., << 1, so that ub and the magnetic pressure pb can be practically neglected in 
the fiow equations.) Since the nonthermal particles are relativistic their pressure is given 
by Pe = Mg/S. The pressure of the thermal plasma is related to its energy density through 
Pgas = (7 — ^)ugas, whcrc 7 is the adiabatic index associated with the thermal gas that we 
take to be either 7 = 4/3 or 7 = 5/3, depending upon whether the temperature of the 
thermal plasma is relativistic or not. The enthalpy of the shocked gas is given by, 

Wjs = Uj^mpC^ + 7/(7 - l)P9as + 4pe- (15) 

The total pressure of the shocked layer is the sum pjg = Pe +Pgas, and must satisfy Pjs{,z) = 
Pext{z) by virtue of momentum balance across the contact discontinuity surface (BL07). With 
the above parametrization for u^. we obtain, using Eq. 0151) . 

37^e 

Ueiz) = — 77 rrPextiz). (16) 

Eqs. ( [TT]) - (fT6l) augmented by the continuity equation (eq. 20 in BL07), the normal compo- 
nent of the momentum equation (eq. 26 in BLOT), and the condition Pjs{z) = Pext{z) form 
a closed set for the unknown variables of the shocked fluid. 
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